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Abstract 

Within a thermal model we estimate possible multiplicities of scalar glueballs in central Au+Au 
collisions at AGS, SPS, RHIC and LHC energies. For the glueball mass in the region 1.5-1.7 GeV, 
the model predicts on average (per event) 0.5-1.5 glueballs at RHIC and 1.5-4 glueballs at LHC 
energies. Possible enhancement mechanisms are discussed. 
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Glueballs are probably the most unusual particles predicted by the QCD but not found 
experimentally yet. It is believed that they can be produced in "gluon-rich" processes, 
like the J/ip radiative decay or pp collisions The lowest glueball states with quantum 
numbers J'^'" = 0++, 2++, 0^+ predicted by the lattice calculations lie in the range 

1.4-2.4 GeV 

m(0++) = 1.6±0.2GeV, m (2++) = 2.2 ± 0.1 GeV, m(0"+) = 2.4 ± 0.05 GeV . (1) 

The uncertainty in masses is mainly introduced by the mixing with the qq states. For 
example, it is expected that the /o(1500) meson has a very large glueball component (see 



recent discussion in Ref. j4|). We think that relativistic heavy-ion collisions might be a 
right tool to produce and study these exotic particles. If the quark-gluon plasma is indeed 
produced in such collisions, a significant fraction, more than 30%, of its entropy should be 
represented by thermal gluons. If these gluons survive until the hadronization stage, they 
will inevitably form glueballs as two-gluon bound states. 

It is well established now jsHo,!?! that ratios of hadron multiplicities, observed in central 
heavy-ion collisions in a broad range of bombarding energies, can be well reproduced within 
a simple thermal model. This model assumes that all hadrons are formed in a common 
equilibrated system characterized by the temperature T, the chemical potential /i^ and the 
volume V . The number density N/V of bosons in a baryon-free system is given by the 
formula {% = c = 1) 



V (27r)3 I 



(2) 



exp 

where m and J are, respectively, the boson mass and spin. Following Refs. P, 0] we take 
into account possible deviations from chemical equilibrium for hadrons, containing nonzero 
number (us) of strange quarks and antiquarks, by introducing a strangeness suppression 
factor 

In our previous paper we have used this model to predict possible yields of exotic 
baryonia. In this letter we apply the same model with the same parameters to estimate 
multiplicities of glueballs, assuming that they are in thermal and chemical equilibrium with 
other hadrons. In order to eliminate unknown volume we consider the ratio of the glueball 
multiplicity to the multiplicity of 0(1020) mesons (n^ = 2). Assuming further that the 
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glueball is a flavor-neutral particle we get for this ratio: 



2 J 



■7, 



-2 



3/2 



(m^ - mG)/T 



(3) 
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where mc (^</>) is the glueball ((/)-meson) mass. The r.h.s. of Eq. (jH)) is obtained in the 
Boltzmann approximation and in the lowest order in T/m. Below we consider only the 
scalar glueballs (J = 0). 

TABLE I: Parameters of thermal model for central Au+Au and Pb+Pb collisions at different cm. 
energies and the observed multiplicities of i;^-mesons. 



reaction 



Au+Au 
Pb+Pb 
Pb+Pb 
Pb+Pb 
Au+Au 
Au+Au 
Au+Au 



y/s^ (GeV) 



4.87 
8.87 
12.4 
17.3 
130 
200 
5.5 • 10^ 



r(MeV) 



119.1 5| 

145.5 SI 

151.9 5| 
154.8 

n 

176 6] 

177 6] 
177 



Is 



0.763 5 1 
0.807 SI 
0.766 [51 
0.938 
1.0 6| 

1.0 bi 

1.0 




38+3 



We use the parameters of thermal model obtained in Refs. P| (flt B) and by fltting 
the hadron ratios observed in central Au+Au and Pb+Pb collisions at various energies. 
These parameters as well as experimental values of the 0-meson multiplicities in these 
reactions are given in Table HJ In the case of Au+Au collisions at the AGS and RHIC 
energies the total multiplicities are not yet available. For these reactions we use the 
(p/n ratios observed at midrapiditv \\\ 12, [3| and estimate A^^ multiplying these ratios 
by charged pion multiplicities lisjl extrapolated to the whole rapidity space. At the 
LHC energy y^J/vw = 5.5 TeV we take the same T and 7^ as for central Au+Au collisions 
at ^snn = 200 GeV and use the value A^^ = 100, obtained by extrapolating the RHIC 
data. In this extrapolation we assume that the ratio of multiplicities at the LHC and 
RHIC energies equals the corresponding ratio of charged particle multiplicities, Nrh. The 
phenomenological formula for the energy dependence of Nch is taken from Ref. [l6[. 

By using Eq. Q with the parameters from Table H] we calculate Nq/N^ for the case of 
most central Au+Au collisions at different bombarding energies. The results for two values 




FIG. 1: Ratio of mean multiplicities of glueballs and (/)-mesons in central Au+Au collisions at 
different cm. bombarding energies calculated within the thermal model. Solid and dashed lines 
show results for glueballs with masses mc = 1.5 and 1.7 GeV, respectively. 

of the glueball mass, mc = 1.5 and 1.7 GeV, are shown in Fig. ^ A nontrivial behavior of 
the Nc/N^-TSitio is explained by the nonmonotonic energy dependence of 7^. One can see 
that for the most central Au+Au collisions at ^/sNN ~ 200 GeV this ratio is predicted in 
the range of 1.5-4%. 

Figure |21 shows the excitation functions of the absolute glueball multiplicity Nq in the 
central Au+Au collisions. At higher RHIC energy the predicted glueball yields are in the 
range of 0.6-1.5 (per event) for rriG =1.5-1.7 GeV. At the future LHC facility we predict, 
respectively, 1.5-4 glueballs per central Au+Au collision. 

It is instructive to compare these predictions with glueball multiplicities obtained for 
superposition of independent nucleon-nucleon collisions. In order to estimate glueball yields 
in a single NN-collision we use the thermal model of pp and pp reactions, suggested in 
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FIG. 2: Mean multiplicities of glueballs in central Au+Au collisions for different values of mass mo ■ 

Ref. Within this model, the hadron multiphcities observed in pp interactions at y/s = 
200 GeV can be reproduced with the effective temperature T ~ 175 ± llMeV and the 
particle emitting volume ^ ~ 35 ± 14 fm^ . Using Eq. (j21), we get the following estimate for 
average multiplicity of glueballs in a single pp-collision at ^/s = 200 GeV: 

Ng (pp) ~ (3.4 - 9.0) ■ 10"^ . (4) 

The two values in the r.h.s. correspond to = 1.5 and 1.7 GeV. We further assume 
that probabilities of the glueball production in the pp- and the NN-collisions are the same 
at high energies. Multiplying Eq. (0} by A^part/2 ^178 [18], the average number of NN- 
pairs participating in most central Au+Au collisions, we obtain 0.6-1.6 glueballs per central 
Au+Au collision at ^/snn = 200 GeV. From Fig. |21one can see that these values are rather 
close to the glueball multiplicities, predicted by the thermal model for the highest RHIC 
energy. 

We think, that the estimates presented above should be regarded as a lower bound for 



the glueball yield in high-energy heavy-ion coUisions. As follows from the QCD-based 



calculations [19|, |20|| and also from some effective models |21|, |23, the gluon-like excita- 
tions acquire a temperature-dependent mass, nig ~ 0.8 GeV just before the hadronization 
transition. Such quasiparticles may form glueball-like bound states already in the plasma 
phase j^]. Applying the same formula © for gluonic quasiparticles and replacing 2 J + 1 
by the effective statistical weight of massive gluons, Ug ~ 10, one gets the ratio of the scalar 
glueball multiplicity to the gluon multiplicity at T = Tc. 

3/2 



Ng ^ \mg) 



(5) 



where the numerical value is obtained assuming Tc = 170 MeV and mc = 2mg . 

However, in a rapidly expanding system the abundances of different species will be de- 
termined not only by the temperature, but also by the corresponding reaction rates. In this 
case the Nc/Ng-idXio may significantly overshoot the "equilibrium" value (IS))j2^. Indeed, 
during the hadronization, the massive gluons may decay into gg-pairs or recombine into the 
gg-hormdi. states forming, respectively, mesons or glueballs. Using simple kinetic equations 
for glueball and gluon abundances with the recombination term, proportional to the gg G 
cross section, agg, and the g ^ qq decay term, characterized by the width Tg , we obtain the 
estimate 

Ng ngaggVj-ei 



N r 



70 . 



(6) 



Here we have used the typical values agg ~ 10 mb, Vg ~ 100 MeV, v^^i ~ yTc/fng ^ 0.5 and 
the equilibrium gluon density ng = ng {Tc) ^ 0.05 fm~^. The ratio (jHl) exceeds the thermal es- 
timate (jni) by more than one order of magnitude. Therefore, much more favorable conditions 



24|. 



for the glueball production may be realized in the explosive hadronization scenario 

For the experimental identification of glueballs one can use the decay mode G — >■ KK 
with predicted partial width in the range 10-40 MeV. For typical conditions, when glueballs 
are slow in the cm. frame, one should look for the back-to-back correlated kaons with 
individual momenta of about 0.6-0.7 GeV/c. The G — > 77 channel with the partial width 
of about 1-10 keV can also be used for the glueball identification. 

In conclusion, we have used a simple thermal model to predict possible yields of scalar 
glueballs in central Au+Au collisions at different energies, from AGS to LHC. Their maximal 
multiplicities are predicted in the range of 1.5-4% of the 0-meson multiplicity. Even larger 
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yields are expected in the case of explosive hadronization of the quark-gluon plasma. We 
believe that such yields are in the reach of experimental observations. 
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